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Abstract 

Isotopic yields for light particles and intermediate mass fragments have been 
measured for ii2Sn+ii2Sn, ii2sn+i24sn, I24g^^ii2s^ ^^^^ 1243^^1243^ 

tral collisions at E/A=50 MeV and compared with predictions of stochastic 
mean field calculations. These calculations predict a sensitivity of the isotopic 
distributions to the density dependence of the asymmetry term of the nuclear 
equation of state. However, the secondary decay of the excited fragments 
modifies significantly the primary isotopic distributions and these modifica- 
tions are rather sensitive to theoretical uncertainties in the excitation energies 
of the hot fragments. The predicted final isotope distributions are narrower 
than the experimental data and the sensitivity of the predicted yields to the 
density dependence of the asymmetry term is reduced. 
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Introduction 

The density dependence of the asymmetry term in the nuclear equation of state (EOS) 
is an important but poorly constrained property of nuclear matter [1-3] . Nuclear structure 
data provide few constraints on the density dependence of the asymmetry term [3]. On 
the other hand, the asymmetry term and its density dependence govern the density, radius 
and proton fraction of neutron stars [2] and provide strong motivations for theoretical and 
experimental investigations of these issues. Recently, a number of calculations have tried to 
identify experimental observables, which can provide constraints on the density dependence 
of the asymmetry term [4-7]. 

In this paper, we focus on sensitivities that have been predicted for observables in ener- 
getic central nucleus coUisions [4-7] . At incident energies greater than about 35A MeV, the 
central density in such collisions initially increases as the projectile and target nuclei over- 
lap and then decreases as the system collectively expands and eventually multi-fragments. 
Previous investigations have shown that excited systems produced in such collisions undergo 
bulk multifragmentation characterized by a short breakup time scale 100 fm/c [8-11] and 
final states containing more than four fragments of charge Z > 2 [12,13]. 

For heavy systems in which the neutron density exceeds the proton density, the asymme- 
try term is repulsive for neutrons and attractive for protons. The asymmetry term therefore 
enhances the dynamical emission of neutrons relative to protons in such collisions; the degree 
of enhancement reflects the magnitude of the asymmetry term and its density dependence 
[4-7]. The difference between neutron and proton emission rates in such collisions can either 
be probed by direct measurements of pre-equilibirum neutron and proton spectra or by ex- 
amining the isotopic composition of the bound fragments that remain after emission [4-7]. 
In this paper, we will concentrate on the fragment observables. 

Fragment observables for these reactions have been described successfully via either sta- 
tistical [14-17] or dynamical [18,19] models. To investigate the dependence of such observ- 
ables on the density dependence of the asymmetry term, it is necessary to calculate the 
relative emission of neutrons and protons and assess the change in the isotopic composi- 
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tion in the prefragment during the expansion stage prior to the multifragment breakup [4] . 
Within the context of a hybrid model, this pre-equihbrium emission was calculated by ref. 
[20] using a Boltzmann-Uehling-Uhlenbeck (BUU) equation model of ref. [4] and the sub- 
sequent fragmentation was explored using the equilibrium Statistical Multifragmentation 
Model (SMM) of ref. [21]. In a second approach [22], the expansion and multifragmentation 
of the system was calculated by a statistical rate equation approach for surface emission 
called the Expanding Emitting Source (EES) model of ref. [14] . Both of these calculations 
predicted that the final isotopic composition of observed fragments should be sensitive to 
the density dependence of the asymmetry term. 

The predicted sensitivities for the surface emission in the EES calculations and for 
the bulk emission in the hybrid BUU-SMM calculations, however, are completely oppo- 
site [20,22]. For calculations using an asymmetry term with softer density dependence, the 
EES approach predicts preferential surface emission of more neutron rich fragments while 
the BUU-SMM approach predicts preferential bulk emission of more isospin symmetric frag- 
ments. These differences stem from different model assumptions in these two approaches 
about the density distribution of the system at the time of fragment production. In the EES 
approach, the fragments are emitted at normal density along with the protons and neutrons 
from a residue, which is at subnuclear density. In the BUU-SMM approach, the fragments 
originate from the bulk disintegration of the residue itself. 

In this paper, we investigate how the isospin transport and dynamics is related to the 
asymmetry term within the dynamical stochastic mean field theory (SMF) approach wherein 
the evolution of the density and nuclear mean field is calculated self consistently [7,23]. In 
previous publications, stochastic mean field theory predictions for the infiuence of the asym- 
metry term on fragment production, collective fiow, incomplete fusion and binary collisions 
have been reported [24-27]. Here, we compare this model to isotopically resolved multifrag- 
mentation data measured in central ■'^^^Sn-|--'^^^Sn and ■'^^^Sn-|--'^^^Sn reactions at 50A MeV. 
These experimental multifragmentation data are presented in the next section. This is fol- 
lowed by a description of the SMF approach, which provides predictions for the dynamical 
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production of highly excited fragments, and of the decay of these excited fragments via the 
MSU statistical decay code [28,29]. The paper concludes with a discussion of the comparison 
between data and theory and the issues that remain for future investigations. 

The experiment 

Central ^^'^Sn+^^'^Sn, ii2Sn+i24g^^ i24gj^^ii2gj^ i24g^^i24gj^ collisions were measured 
at the National Superconducting Cyclotron Laboratory at Michigan State University, using 
5 mg/cm^ ^^^Sn and ^^^Sn targets and 50 MeV per nucleon ^^^Sn and ^^^Sn beams. Isotopi- 
cally resolved light particles and intermediate mass fragments with 3<Z<8 were measured 
with the Large Area Silicon Strip detector Array (LASSA) [30,31], an array consisting of 
9 telescopes, each comprised of one 65 /im and one 500 /im Si strip detector, followed by 
four 6-cm thick CsI(Tl) detectors. The 50mm x 50mm lateral dimensions of each LASSA 
telescope are divided by the strips of the second silicon detector into 256 (3x3 mm2) square 
pixels, providing an angular resolution of about 0.43''. The LASSA device was centered at 
a polar angle of Oiab—^'^" with respect to the beam axis, providing coverage at polar angles 
of 7° < 6iab <58°. At other angles, charged particles were detected in 188 plastic scintillator 
- CsI(Tl) phoswich detectors of the Miniball/ Mini wall array [32], which subtended polar 
angles 7° < Oiab <160°. The Miniball/ Miniwall array provided isotopic resolution for H and 
He nuclei and elemental resolution for intermediate mass fragments (IMF's) with 3<Z<20. 

The total charged particle multiplicity detected in the two arrays was used for impact 
parameter determination. Central collisions, corresponding to a reduced impact parameter 
of b/bmax <0.2 [33], were selected by a gate on the top 4% of the charged-particle multi- 
plicity distribution. From cross section measurements for such events, we estimate a value 
for bmax— 7-51 fm by neglecting fluctuations and assuming that multiplicity decreases mono- 
tonically with impact parameter. This would suggest that our impact parameter selection 
corresponds to b <1.5 fm; however, multiplicity fluctuations at fixed impact parameter may 
extend the included range of impact parameters outward (up to b ~3 fm). 

In the following we present isotopically resolved differential multiplicities for fragments 
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emitted at center of mass angles of 70" < 6cm <110°. At these angles, the coverage of 
the LASS A array is excellent; the only losses occurred for fragments emitted at very low- 
energies E/A< 0.2 MeV in the center of mass, corresponding to small laboratory angles of 
Giab ~7°. This enabled accurate calculation of the detection efficiency for 70" < 9cm <110''; 
the fragment spectra were fitted and the losses at low energies (<2%) were estimated and 
corrected. The data presented below have been corrected for the losses below threshold, for 
inefficiencies in the solid angle coverage and for multiple hits in the detector telescopes. 

Figure 1 shows the measured average differential multiplicities of Li, Be, B, C, N and 
O isotopes at 70° < 9cm <110°. In this figure, the solid squares and circles show the 
ii2Sn+ii2Sn and i24Sn+i24Sn data, respectively. The isotopic yields of ii2Sn+i24Sn and 
^^^Sn+^^^Sn are essentially equal; they have been averaged and are shown by the open 
diamonds. The x-axis, N-Z, corresponds to the neutron excess of the nuclides. The peaks 
of the distributions are located at isotopes with N=Z+1. The yields of the B, C isotopes 
are offset by a factor of 10 and the yields of the N and O isotopes are offset by a factor of 
100 in the figure. As expected, more neutron rich nuclides are produced by the neutron rich 
system, ^^^Sn+^^^Sn, while the opposite is true for emission from the proton-rich isotope 
system, -"^^^Sn+^^^Sn. The experimental results indicate that the multiplicities of IMF's are 
approximately 10-20% larger for the ^^^Sn-|-^^^Sn entrance channel than for the ^^^Sn-|-^^^Sn 
entrance channel, consistent with previous observations at an incident energy of 40 MeV per 
nucleon [34]. 

In general, the drop from the peak toward more proton-rich isotopes is rather steep 
especially for elements with even values of Z. The main differences between the isotope 
yields for the four different systems are observed in the tails of the isotope distributions, 
where it is greater than a factor of 4 for ^°0. Larger differences may be expected for even 
more exotic isotopes, but the background in the present measurement due to multiple hits 
in the LASSA telescopes does not allow for their accurate determination. 

Recently, it has been shown [35] that the isotopic yields for systems produced at approxi- 
mately the same excitation energy per nucleon or the same temperature satisfy an isoscaling 

6 



relationship. Specifically, the ratio R2i(N, Z) = Y2(N,Z)/Yi(N, Z) constructed using the 
isotope yields Yj(N, Z) with neutron number N and proton number Z from two different 
reactions denoted by the index i, i=l-2, obeys a simple relationship [20,22,35,36]: 

R2i{N,Z)= Ce"^+'^^ (1) 

Here, C is an overall normalization factor and and are isoscaling parameters. This parame- 
terization is discussed in greater detail within the isoscaling section below. If we adopt the 
convention that reaction 2 is more neutron- rich than reaction 1, one expects a to be positive 
and (3 to be negative. We have adopted that convention here and have fitted the ratios of 
the isotopic yields for these four systems to extract the corresponding values for a and /3. 
These values for and are given in Table 1. 

Equation 1, with only three parameters, C, a and f3 can be used to predict the isotope 
yields of ^^^Sn-F^^^Sn as well as the mixed systems, ^^^Sn-F^^'^Sn or ^^^Sn-F^^^Sn, using the 
measured yields of one system, e.g. ^^^Sn-|-^^^Sn. To illustrate how well this parameterization 
relates the yields of these four systems, we take the yields of the ^^''Sn-l-^^^Sn system as a 
reference and use those yields and the fitted values of a and /3 to predict the yields for 
the other four systems. The dash and dot-dashed hues in Fig. 1 are the calculated yields 
for ^^^Sn-|-^^^Sn and ^^^Sn-|-^^^Sn, respectively. For this limited range of asymmetry, these 
isoscaling parameters can be described by a linear dependence on either the initial N/Z or the 
asymmetry parameter, 5 — of the reactions [35]. The excellent agreement between the 
predicted yields and the data suggests that such scaling law extrapolations may have useful 
predictive power. For example, we expect that these scaling predictions can be accurately 
extrapolated to other mass-symmetric systems of A=200-250 nucleons at the same incident 
energy per nucleon but with very different isospin asymmetry. 

Description of the Models 

Now we turn to the theoretical interpretation of these data. To study the density de- 
pendence of the asymmetry term of the EOS, we adopt the viewpoint of the Stochastic 
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Mean Field (SMF) approach described in refs. [7,23]. In this approach, the time evolution 
of the nuclear density is calculated by taking into account both the average phase-space 
trajectory predicted by the Boltzmann-Nordheim-Vlasov equation and the fluctuations of 
the individual collision trajectories about this average that can be predicted by equations 
of the Boltzmann-Langevin type. 

The virtue of such a dynamical approach for the study of isotopic effects lies in its self- 
consistency. The flow of neutrons and protons is calculated under the influence of Coulomb 
and asymmetry terms, which reflect self-consistent ly the motion of these nucleons. Several 
different density dependences of the asymmetry term were explored from which two are 
selected for presentation here. In both cases, the asymmetry term is approximated by the 
form, 

Esymip,S)=S{p)-S^ (2) 

where for the asymmetry term with a stronger density dependence, 

S(p)^a.(P-Y\b. ^^PlP^}\ . (3) 
\PoJ l + {p/po) ^' 

Here, p is the physical density, po the saturation density, a=13.4 MeV and b=19 MeV 
[2,4,37]. In the following, we refer to this as the "super stiff" asymmetry term. For the 
asymmetry term with weaker density dependence, 

S{p)^a-iyj + 240.9p -819.1p2 (4) 

where a=12.7 MeV [25], we refer to this as the "soft" asymmetry term. In Fig. 2, it can be 
seen that the two expressions are nearly equal at saturation density but differ at densities 
that are either much larger or smaller than pq. 

In addition to the asymmetry term, the SMF calculations have a Skyrme type isoscalar 
mean fleld with a soft equation of state for symmetric matter characterized by an incom- 
pressibility constant K=201 MeV. The isoscalar mean field and the asymmetry term of 
these equations of state are used for the construction of the initial ground state and for the 
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time evolution of the collision. The nucleon-nucleon collisions by the residual interaction 
are calculated from an energy and angle dependent parameterization of the free nucleon- 
nucleon interaction and the isospin dependence of the Pauli-blo eking is considered during 
these coUisions. 

The calculation solves the transport equations by evolving test particles of finite width. 
As mentioned above, we use a reduced number of test particles (50 test particles per nucleon) 
in the present calculations to inject numerical noise into the evolution. In test calculations, 
we alternatively employed the fluctuation mechanism discussed in ref. [23], which involves 
damping the numerical noise by utilizing a large number of test particles and introducing 
explicitly physical noise according to thermal fluctuations. It was checked that both methods 
lead to similar results. In contrast to the BOB method of ref. [38], these methods of inserting 
fluctuations are well suited to reactions at flnite impact parameter because they do not 
presuppose knowledge of the most unstable modes. 

When the system expands and reaches the spinodal instability (after about 110-120 
fm/c), the most unstable modes are amplified and initiate the formation of fragments via 
spinodal decomposition. The evolution of the system is continued after spinodal decom- 
position until freeze-out where the number of dynamically produced fragments and their 
properties are finally determined. The system is decomposed into fragments using essen- 
tially a coalescence mechanism in coordinate space; specifically, fragments are defined by 
regions of density in the final distributions that are above a "cut-off" density of l/8po- By 
definition, the freeze-out time occurs when the average calculated number of fragments sat- 
urates. This occurs about 260 fm/c after initial contact of projectile and target nuclei in the 
present simulations. The excitation energy of the fragments is determined by calculating the 
thermal excitation energy in a local density approximation. The procedure is rather rough 
and will overestimate the excitation energy particularly for light fragments. 

Some of the important features of these calculations and of the prior BUU-SMM and 
EES calculations can be understood simply by considering the infiuence of the density de- 
pendence of the asymmetry term on the relative emission rates of neutrons and protons. In 
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all Sn+Sn collisions, the symmetry energy in the liquid drop model is positive, i.e. repulsive. 
The interaction contribution to the symmetry energy gives rise to a repulsive contribution 
to the mean field potential for neutrons and an attractive contribution to the mean field 
potential for protons. The mean field potential for an asymmetry term with stronger density 
dependence is larger at high density and smaller at low density than that for an asymmetry 
term with weaker density dependence. 

It is the low-density behavior that dominates the predictions for the isoscaling param- 
eter. As the system expands and eventually multi-fragments, the prefragment remains at 
subnuclear densities for a long time while it is emitting nucleons. The asymmetry term 
with weaker density dependence around po increases the difference between the neutron 
and proton emission rates leading to a more symmetric prefragment than is produced by 
calculations with the asymmetry term which has a stronger density dependence. 

The SMF calculations are interesting because they are free, in principle, of arbitrary 
assumptions about whether the fragments are formed at the surfaces or from the bulk dis- 
integration of the system. Comparisons between the fragmentation dynamics for different 
asymmetry terms were reported in ref. [39]. The trends of these calculations are consis- 
tent with the prefragment isospin dependences discussed above. In particular, fragments 
produced in calculations with an asymmetry term with strong density dependence tend to 
be more neutron-rich than the fragments produced in calculations with an asymmetry term 
with weak density dependence. In this respect, these predictions are similar to the results of 
the BUU-SMM calculations of ref. [20] and opposite to the results of the EES calculations 
of ref. [22]. 

The SMF fully dynamical formation of fragments should actually be more sensitive than 
the hybrid BUU-SMM to the interplay of the EOS, i.e. to the density dependence of the 
asymmetry term, with the fragmentation process. In the hybrid BUU-SMM calculations, 
the EOS is entering only in the "pre-equilibrium" nucleon emission described above. In the 
SMF approach, wc have not only this isospin effect on fast particle emission but also the full 
dynamics of the isospin fractionation mechanism during the cluster formation. In a neutron- 
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rich system, this leads to a different N/Z "concentration" in the liquid phase (the fragments 
are more symmetric) and in the gas phase (nucleons and light ions, are more neutron-rich) 
[16,40]. This effect is associated with the unstable behavior of dilute asymmetric nuclear 
matter and so in this way we have the chance of testing the EOS also at sub-saturation 
density. 

Asymmetry terms with weaker density dependence around po must show a faster increase 
at low densities and so a larger isospin fractionation/distillation during the fragment forma- 
tion [7] . Therefore in a fully dynamical picture of fragmentation events a "soft " behavior of 
the asymmetry term around saturation density will enhance the formation of more symmet- 
ric fragments for two converging reasons: i) A larger " preequilibrium" neutron emission rate 
as discussed before; ii) A stronger isospin fractionation/distillation during the bulk disinte- 
gration. Opposite effects are of course predicted for a rapidly increasing ("stiff") asymmetry 
term around po- In this sense we can expect the SMF results to be more sensitive to the 
isospin dependences of the EOS at sub-saturation density. 

At freezeout, the fragments are highly excited. For simplicity, we assume that the de- 
excitation of these fragments can be calculated as if the fragments are isolated. For this de- 
excitation stage, we have tabulated the data table states, spins, isospins and branching ratios 
for nuclei with Z<15. Where experimental information is complete, it is used. Alternatively, 
empirical level density expressions are used for the discrete levels. These discrete levels are 
matched to continuum level density expressions as described in ref. [28]. The decay of 
primary fragments with Z>15 are calculated, following ref. [28], using known branching 
ratios, when available, and using the Hauser-Feshbach formalism when the information is 
lacking. The decays of heavier nuclei are calculated using the Gemini statistical decay code 
[41]. 

While the SMF calculations predict the numbers and properties of the hot fragments that 
are produced at breakup, the predictions for the relative abundances of light clusters such as 
the isotopes with Z=l-2 that are emitted before the system expands to sub-nuclear density 
are not very realistic. This prevents a precise modeling of complete events including their 
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detection efficiency and means tliat tfie impact parameter selection based on multiplicity 
cannot be imposed straightforwardly on the calculated events as on the data. This and 
the considerable numerical effort it requires have persuaded us to limit our comparisons to 
calculations composed of 600 events for each of the ^^^Sn+^^^Sn and ^^^Sn+^^^Sn reactions at 
a fixed impact parameter of h—2 fm. We note, however, that the widths in the multiplicity 
distributions at fixed impact parameter are large enough that a range of impact parameters 
may contribute significantly to the experimental data. Future calculations will be necessary 
to assess quantitatively the importance of this impact parameter smearing. 

Overall behavior predicted by the SMF calculations 

In Fig. 3, the solid circles and open squares in the left panel show the measured elemental 
multiphcities for 2<Z<8 averaged over 70" < Ocm <110° for ^^^Sn+i^^gn and ^^^Sn 
collisions, respectively. The right panels show the corresponding measured multiplicities as 
a function of the fragment mass. These averaged multiplicities were obtained by summing 
the isotopic multiphcities for 2<Z<8. 

The dashed lines denote the corresponding distributions of hot primary fragments calcu- 
lated by the SMF model using the super stiff EOS. Due to the low total number of events, we 
averaged these calculations over a slightly larger angular interval of 60° < 9cm <120°. The 
solid hues show the multiplicities of cold fragments after secondary decay. The statistical 
uncertainties in these calculations are shown in the figure as vertical bars. The correspond- 
ing uncertainties in the data are smaller than the data points. If the angular integration 
was performed over the entire solid angle, the averaged calculated multiplicities are about 
20-30% larger. This difference reflects an anisotropy in the calculated primary angular dis- 
tributions for the heavier fragments. In the present calculations, however, we do not have 
the capability to accurately calculate the modifications of the angular distribution due to 
secondary decay so we presently cannot explore this issue more quantitatively. As we will 
show later, this anisotropy has no impact on the shapes of the isotopic distributions for 
Z=3-8. 
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In general, the calculated primary and secondary fragment multiplicities are smaller than 
the measured values for the lighter fragments Z—3,4 and are somewhat closer to the mea- 
sured values for Z=6-8. The lighter fragments with Z<4 are mainly produced in secondary 
decay stage of the theoretical calculations; the primary yields of these hght fragments are 
much smaller relative to the final yields than are the values for the heavier fragments. Be- 
cause the fragment multiplicities and angular distributions depend on impact parameter, 
the comparison shown in Fig. 3 may be sensitive to the impact parameter ranges included 
in both calculation and data. Future calculations over a wide range of impact parameters 
are needed to address this issue. Concerning the greater discrepancy for Z=2-4 fragments, 
we have already noted that the formation of light clusters in the dynamical stage before 
breakup is not well described in BUU- and SMF-type simulations, because the unique struc- 
tural properties of these fragments are not therein well treated. (Treatments of the emission 
of light clusters in coupled transport equations for nucleons and light clusters can be found 
in refs. [42,43] and in the framework of FMD [44] or AMD [45] simulations.) On the other 
hand, there is a considerable emission of protons and neutrons during this stage; the total 
emission and consequently the asymmetry of the remaining source may still be realistic. 

Now we turn to an examination of calculated isotopic yields. The upper left panel of 
Fig. 4 shows the isotopes of carbon nuclei predicted by the SMF calculations over the entire 
angular range for ^^'^Sn-l-^^'^Sn (solid line) and ^^^Sn-|-^^^Sn (dashed line); the dotted-dashed 
and dotted lines show the corresponding calculations over the 60" < 9cm <120° gate. Not 
surprisingly, the more neutron rich ^^^Sn+^^^Sn system preferentially produces the more 
neutron rich isotopes. The peak of the carbon primary distribution for the ^^^Sn-|-^^^Sn 
system occurs at about ^^C while the peak for ^-'^^Sn-l-^-'^^Sn system occurs at lower mass, i.e. 
somewhere between ^^C and ^"^C The differences between the angle gated and total primary 
yields are small, and these difference translate into negligible differences in the shape of 
the isotopic yield distribution after secondary decay; we therefore do not plot the gated 
data because the two curves are indistinguishable when normalized to each other. As the 
statistics make it difficult to perform comparisons to isotopic yields with a 60° < 9cm <120° 
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gate imposed, the remaining calculated multiplicities in the paper are integrated over the 
entire solid angle. 

After sequential decays, one obtains the secondary distributions shown in the lower left 
panel. No longer is there a noticeable difference between the peak locations (at ^^C in both 
systems); instead, the main differences are found in the shape of the distribution, which 
is higher in the neutron- rich isotopes and lower in the neutron-deficient isotopes for the 
^^'^Sn-l-^^^Sn system than it is for the ^^^Sn-|-^^^Sn system. Such trends are also quahtatively 
observed in the experimental data shown for the ■^^^Sn-|--'^^^Sn system by the sohd circles 
and for the ^^^Sn-|-^^^Sn system by the open squares in the lower left panel. However, the 
experimental distributions are considerably wider and more neutron rich than the model 
predictions. This trend is rephcated in the elemental distributions for all of the other 
measured elements. 

Another way to quantify the differences in the isotope distributions is by the asymmetry 
parameter 6—{N-Z)/(N+Z). The average asymmetry < 5 > of the isotopic distribution for 
each element is shown as a function of Z in the right panels of Fig. 4. Following the same 
convention as in the left panels, the solid and dashed lines show the average asymmetries 
for ^^^Sn-|-^^^Sn and ^^^Sn-|-^^^Sn collisions; the upper and lower panels present results for 
the primary and secondary fragment distributions, respectively. The calculated differences 
between the two systems are more pronounced prior to secondary decay than afterwards. 
The asymmetries of the corresponding data, shown for the ^^^Sn-|-^^^Sn system by the solid 
circles and for the ^^^Sn-|-^^^Sn system by the open squares in the lower right panel, are 
larger and display a stronger dependence on the asymmetry of the system than do the final 
calculated fragment yields after secondary decay. 

Isoscaling analyses 

A more sensitive way to compare isotopic distributions is to construct the isotopic ratio 
R2i(N,Z)=Y2(N,Z)/Yi(N,Z) from the isotope yields Yj(N,Z) with neutron number N and 
proton number Z from two different reactions. As discussed in the experimental section. 
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R2i(N,Z) obeys a simple relationship R2i{N,Z) = C e'^^+^^ where C is an overall nor- 
malization factor and a and /3 are isoscaling parameters [20,22,35,36]. Such an isoscaling 
relationship can be obtained in statistical theories for two systems that are at the same 
temperature when they produce fragments. Binding energy factors common to the yields 
for the fragments in each system are cancelled by the ratio when the temperatures are equal, 
leaving terms related to the chemical potentials or separation energies [22] . In grand canon- 
ical models of multifragmentation a — and (3 — for example, where n,p are the 
differences in the chemical potentials for the neutrons and protons in the two systems and 
T is the temperature [20,22,35,36]. In some calculations [36], the values for the isoscaling 
parameters extracted from equilibrium multifragmentation models are similar before and 
after sequential decays, an observation that has been attributed to a partial cancellation of 
secondary decay effects [35, 36]. 

While isoscaling can be expected for many statistical processes [20,22,35,36], the ques- 
tion of whether it can be expected for specific dynamical calculations remains open. To 
investigate whether the SMF dynamical model displays isoscaling, we construct the relative 
isotope ratios, R21, primary, using the primary fragments produced in ^^^Sn-|-^^^Sn collisions 
as reaction 2 (numerator) and in ^^^Sn-|-^^^Sn collisions as reaction 1 (denominator). 

The results are shown in the upper panel of Fig. 5. The error bars reflect the statistical 
uncertainties. The predicted isotope ratios for these primary fragments depend very strongly 
on the neutron number and follow trends that appear consistent with isoscaling relationship 
defined by Eq. 1. The uncertainties are large reflecting the low statistics of the simulations, 
but the strong dependence on neutron number makes it possible to discern apparent isoscal- 
ing trends nonetheless. The lines are best fits using Eq. 1 resulting in C=0.96, q;=1.07 
and /9=-1.43. These values for a are much larger than values observed in the experiment. 
The lower panel provides the corresponding SMF predictions for the ratios, ^2ijinai-i of the 
yields of particle stable nuclei after secondary decay. For comparison purposes, the scale for 
the ordinates of the top and bottom panels are chosen to be the same; this demonstrates 
graphically that the trends of the final isotope ratios are much fiatter and the corresponding 
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isoscaling parameters (q;=0.286 and (3— -.288) are much smaller. Clearly, the isoscahng pa- 
rameters predicted by dynamic SMF calculations are strongly modified by secondary data. 
This trend is very different from equilibrium statistical models for multifragmentation where 
the isoscaling parameters have been predicted to be insensitive to secondary decay [35,36]. 

The isoscaling behavior of the dynamically produced fragments arises not from thermal 
physics but rather from some special characteristics of the SMF primary distributions pre- 
dicted for these reactions. We find, for example, the SMF primary isotopic and isotonic 
distributions can be roughly described by Gaussians, see Fig. 3 [46]. Isotopic distributions, 
for example, can be described by 



y(7V,Z) = /(Z)exp 



N -N{Z 



2<jI 



(5) 



where N (Z) is the centroid of the distribution and (7| describes the width of the distribution 
for each for each element of charge Z. This leads to an exponential behavior of the ratio R21, 
since, neglecting quadratic terms in N, 

1 



lni?2i = — 



N{Z),-N{Z), 



N 



(6) 



Note Eq. 6 requires the values for cr| to be approximately the same for both reactions. 
We have observed this to be the case for our SMF calculations of Sn-|-Sn collisions (to 
within the statistical accuracy ~10%). For the ratios for every element, to be optimally 
described by the same parameter, the ratio N {Z)r^ — N (Z)^ /az must be independent of 
Z. The statistics of the calculation do not allow a detailed test of this assumption, but it 
does appear that this ratio increases somewhat with Z, as Fig. 5 confirms. The primary 
distributions therefore do not respect the isoscaling relationship as well the data do. 

At variance with the statistical fragmentation models, the secondary decays substantially 
modify the isoscaling parameter. The width a% decreases due to secondary decay and the 



difference 



N2 — Ni likewise decreases fractionally, but by a larger amount. Moreover, the 



final shape is no longer Gaussian, but due to secondary decay, it refiects the binding energy 
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as a function of neutron excess more strongly (see Fig. 4). These changes combine to 
decrease the isoscahng parameter as shown in the lower panel of Fig. 5. 

Sensitivity of the SMF calculations to the asymmetry term 

The density dependence of the asymmetry term has a significant influence on the relative 
emission rates of the neutrons and protons and, consequently, on the isospin asymmetry of 
the hot fragments prior to secondary decay. As discussed previously, an asymmetry term 
with weaker density dependence tends to remain more important at lower densities, driving 
the fragments closer to isospin symmetry, than does an asymmetry term with stronger 
density dependence. Consistent with this general consideration, the calculated primary 
isotope distributions in ^^^Sn+^^'^Sn collisions, shown in Fig. 6 for carbon (upper left panel) 
and oxygen (upper right panel), are more neutron- rich for the super stiff asymmetry term 
(solid hne) than they are for the soft asymmetry term (dashed line) . A similar trend is also 
predicted for the ^^^Sn+^^^Sn system, but is not shown in the interest of brevity. 

A similar trend is observed in the corresponding final distributions that are obtained 
after secondary decay and shown in the middle panel with the same convention for the solid 
and dashed lines as in the upper panel. Both secondary distributions calculated for super 
stiff and soft asymmetry terms, however, are significantly narrower and more proton-rich 
than the experimental distributions shown by the closed circles in the figure. (The lower 
panels, which display corresponding calculations when the excitation energy is reduced by 
50%, will be discussed in the next section of this paper.) Similar trends are also observed 
for the ^^^Sn-|-^^^Sn and for the other elements with 3<Z<8, though we do not for brevity's 
sake show those results. 

In Fig. 7, we present the related dependence of the SMF predictions for the isotope 
ratios R21 upon the density dependence of the asymmetry term. We take advantage of 
the fact that the results in the Fig. 5 can be compactly displayed by the scaled function 
S (N) — R21 {N, Z) ■ e~^^ which condenses the isotopic dependence for the various elements 
onto a single line [22] . The left panel in Fig. 7 shows the results for the super stiff asymmetry 
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term and the right panel shows the results for the soft asymmetry term. In each panel, the 
values for S(N) obtained from the primary distribution are shown by the symbols clustered 
about the dashed lines, the results obtained from the secondary distribution are shown by 
the symbols clustered about the dot-dashed hues and the results from the data are shown 
by the solid lines in each panel to provide a reference. Both the primary and secondary 
values for S(N) have been fit by exponential functions to obtain corresponding values for 
the scaling parameter and these values are given in the figure. 

Generally, the primary distributions for both equations of state display a much stronger 
dependence on neutron number than do the final isotopic distributions and the data. How- 
ever, the influence on the isoscaling parameter is statistically not very significant. Indeed, 
as we pass from a " stiff" asymmetry term to a " soft" one, we do have a stronger isospin 
fractionation/distillation, as already discussed before. The centroid of the distribution, N2, 
decreases (see Fig. 6) but also the width erf decreases; the decrease in width, however, is of 
the order of 10% and comparable to its statistical uncertainty. The calculated final distribu- 
tions display a weak sensitivity to the density dependence of the asymmetry term; the values 
for (q;=0.286) obtained for the super stiff asymmetry term are larger than the values for 
(q;=0.254) obtained for the soft asymmetry term. The sensitivity to the asymmetry term is 
considerably less than that reported for the EES model [22], and for the BUU-SMM hybrid 
calculations [20] . Unhke these latter two calculations, both super stiff and soft asymmetry 
terms yield a values that are significantly lower than the value extrapolated from the data 
(q;=0.36). One should note, however, that the excitation energies of these latter calculations 
could be more freely varied to achieve better agreement with the experimental observations. 

Discussion and summary 

The calculated final isotopic distributions for both asymmetry terms differ from the 
measured ones in that they are narrower; more neutron deficient; and show a weaker depen- 
dence on the isotopic asymmetry of the total system. The last characteristic is reflected more 
clearly by the isoscaling parameter than by direct examination of the isotopic distributions. 
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themselves. In these respects, the calculated results for the two different asymmetry terms 
are more similar to each other than they are to the data. We believe that it is probably pre- 
mature at this stage to focus attention on the sensitivity of the predicted final distributions 
to the asymmetry term. Instead, let us concentrate upon what may be required to bring 
the final isotopic distributions into greater concordance with the measurements. 

The tendencies of the final isotopic distributions to be more neutron deficient and to dis- 
play a weaker dependence on the isotopic asymmetry of the system are somewhat related. 
Both point to difficulties the present calculations have in producing neutron-rich isotopic 
distributions and indicate a surprising sensitivity of the final results to the primary distri- 
butions. That the final isotopic distributions are too neutron deficient may result from the 
primary distributions being too neutron deficient on the average, too narrow (i.e. u\ is too 
small) or that the secondary decay calculations predict too much neutron emission after 
freezeout because the calculated excitation energies are too high or the excitation energy 
distributions are too narrow resulting in the loss of components at low excitation energies 
that could decay to neutron-rich final products. 

We note that additional experimental measurements may help to resolve these questions. 
The average isospin asymmetry of the initial distributions is trivially related by charge and 
mass conservation to the average isospin asymmetry of the nucleons emitted during the SMF 
calculations before the freezeout (t=260 fm/c) chosen for these calculations. Comphmentary 
measurements of the yields and energy spectra of light particles can help to determine 
whether these missing neutrons are carried away primarily during pre-equilibrium emission 
during the compression-expansion stage or during the later evaporative decay of the hot 
fragments. 

Previous authors have identified issues relevant to our calculations, which may influence 
the asymmetries of the hot fragments at freezeout [43,47,48]. As discussed above, the present 
simulations underestimate the emission of light clusters (d, t, ^He, ^He, ^Li, ^Li, etc.) during 
the dynamical evolution prior to the freezeout. Previous studies [43,47,48] have noted that 
the neglect of the emission of ^He emission is particularly problematic because it is abundant 
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and because each ^He particle enhances the isospin asymmetry of the remaining system by 
removing four nucleons without changing the neutron excess. Indeed, it has been speculated 
that that ^He emission may have an influence on the isospin asymmetry of the other clusters 
and fragments that is of the same order of magnitude as the influence of the mean field 
[43,47,48] . The issue needs additional theoretical attention. 

Concerning neutron emission after freeze-out, we note that the number of neutrons re- 
moved by secondary decay depends primarily on the fragment excitation energies and the 
relative branching ratios for neutron and charged particle emission. There are significant 
uncertainties in the calculation of the excitation energies of the fragments, which are related 
to the difficulty to estabhshing their precise ground state binding energies. To explore the 
sensitivity of the results to the excitation energy, we have reduced the excitation energy of 
each fragment by a multiplicative factor f where 0.5<f<l and recalculated the final fragment 
isotopic distributions. 

The solid and dashed lines in Fig. 6 for carbon fragments (bottom left panel) and oxygen 
fragments (bottom right panel) show the calculated final distributions for f=0.5 using super 
stiff and soft asymmetry terms, respectively. Clearly, it is possible by reducing the excitation 
energy to shift the isotope distribution in the direction of the more neutron-rich isotopes, so 
as to make the mean isospin asymmetry of the calculated final and measured distributions 
to be the same. However, the widths of the calculated final isotopic distributions will still 
be narrower than the measured ones. 

This discrepancy between the theoretical and experimental widths would be reduced if 
the theoretical primary distributions were wider in their excitation energy distributions or 
wider in their isotopic distributions or both. It is interesting that the primary distributions of 
equilibrium statistical model calculations that reproduce the experimental final distributions 
are much wider in excitation energy and neutron number than those predicted by the SMF 
calculations [21,36]. Future investigations will be needed to address whether wider primary 
distributions in excitation energy or neutron number can be attained in the SMF model by 
altering some of the underlying model assumptions such as the manner in which they are 
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defined at freezeout. 

Increased widths may be achieved by performing calculations for a range of impact pa- 
rameters rather than the single impact parameter b=2 fm presented here. We note that 
the inclusion of larger impact parameter events may broaden the primary distributions at 
mid-rapidity because it will require the inclusion of fragments emitted from the neck joining 
projectile- and target-like residues. In ref. [7], it was shown that SMF calculations predict 
such " neck" fragments to be more neutron rich because the isospin fractionation/distillation 
effect is somewhat reduced in peripheral events, leaving more neutrons in the fragments, and 
there can also be an overall neutron enhancement in the neck region for such events due to 
the neutron skins of the projectile and target [7]. Moreover, the excitation energies of the 
primary fragments in more peripheral collisions are also somewhat reduced. Experiments 
also suggest that neck fragments in very-peripheral collisions are more neutron- rich [49] , but 
there is little detailed experimental information about the dependence at smaller impact pa- 
rameters. Combining more calculated peripheral events with the calculated central collision 
events presented in this paper may result in the broader isotopic distributions that appear 
to be required by the data. 

In summary, we have measured the isotope distributions of Z=2-8 particles emitted in 
four different Sn-|-Sn reactions with different isospin asymmetry and have calculated them 
with a dynamical model that includes fluctuations that give rise to fragment production. 
The experimental data display a strong dependence on the isospin asymmetry that can 
be accurately described by an isoscaling parameterization. The theoretical calculations 
reproduce the yields for the heavier fragments with Z=6-8, but underpredict the yields of 
the lighter ones, which are not strongly produced as primary fragments. The calculated 
final isotopic distributions display isoscahng, but the calculated isotopic distributions are 
narrower, more neutron deficient; and show a weaker dependence on the isotopic asymmetry 
of the system than do the data. The density dependence of the asymmetry term of the EOS 
has an effect on the calculated final isotopic distributions. The distributions calculated using 
the asymmetry term with stronger density dependence are more neutron-rich and are closer 

21 



to the measured values. These trends are similar to prior results obtained for a BUU-SMM 
hybrid model, but different from the trends for evaporated fragments predicted by EES 
rate equation calculations. The present level of agreement between theory and experiment 
precludes definitive statements about the density dependence of the asymmetry term of 
the EOS, but it does reveal that that final distributions are surprisingly sensitive to the 
widths predicted by the SMF model for the primary fragment isotope and excitation energy 
distributions. A number of theoretical issues, such as the pre-equilibrium emission of bound 
clusters, the calculations of fragment excitation energies, the way fragments are defined at 
freezeout, and the impact parameter range modeled by the calculation may influence the 
calculated results. Additional theoretical work is required to explore these issues and to 
determine the role they may play in the resolution of these discrepancies. Complementary 
measurements of the isospin asymmetry of light cluster emission prior to the multifragment 
breakup can provide information relevant to the resolution of these issues. 
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TABLES 



R(\u-ti<)ii 2 


Reaction 1 


n 


3 






0.18±0.01 


-0.19±0.01 


i24Sn+i24Sn 


ii2Sn+ii2Sn 


0.36±0.02 


-0.39±0.01 



TABLE L Values for a and (3 for obtained from fitting the isotope ratios R21. 
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FIGURE CAPTIONS: 



Figure 1: Average differential multiplicities at 70° < dcM <110° for Li, Be, B, C, N and O 
isotopes as a function of neutron excess (N-Z) of the isotope. The sohd circles (connected 
by solid lines to guide the eye) are the data for the ^^'^Sn+^^'^Sn system with N/Z=1.48. 
The solid squares are data for the hghtest system ^^^Sn+-'^^^Sn with N/Z=1.24. The open 
diamonds are the averaged values from the two mixed systems, ^^^Sn+^-'^^Sn and ^■'^^Sn+^^^Sn. 
The dashed and dot-dashed lines are predictions from Eq. 1. See text for more details. 
Figure 2: The sohd curve and dashed curves indicate the density dependencies for the super 
stiff and soft asymmetry terms, respectively. 

Figure 3: Differential multiplicities for ^^^Sn+^^^Sn collisions (upper panel) and ^^^Sn+^^^Sn 
collisions (lower panel) as a function of the fragment charge (left panels) and the fragment 
mass (right panels). The points are the data. The dashed and solid lines are the calculated 
primary and final fragment differential multiplicities, respectively. Statistical uncertainties 
are shown for the calculations; the corresponding uncertainties in the data are smaller than 
the data points. 

Figure 4: Left panels: Calculated primary (upper panel), calculated final (lower panel) and 
measured (lower panel) carbon isotopic yields for Sn+Sn coUisions. Right panels: Calculated 
primary (upper panel), calculated final (lower panel) and measured (lower panel) mean 
isospin asymmetries as a function of the fragment charge for Sn+Sn collisions. The lines 
and data points are further explained in the text. 

Figure 5: Upper panel: R21 values obtained from the ratios of the primary isotopic distri- 
butions for ^^^Sn-l-^^^'Sn collisions divided by those for ^^^Sn-|-^^^Sn collisions. Lower panel: 
Corresponding R21 values obtained from the ratios of the final isotopic distributions. Each 
line in the two panels corresponds to ratios for a given element. Elements with Z=2-8 (Z=l- 
8) are represented from left to right in the upper (lower) panel. The lines are the result of 
fitting R21 with Eq. 1; the dependencies on neutron number for the best fits are given in 
each panel. 
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Figure 6: Upper panel: Dependence of the primary distributions for carbon(left panel) and 
oxygen (right panel) upon the density dependence of the asymmetry term. Middle panel: 
Dependence of the final distributions for carbon (left panel) and oxygen (right panel) upon 
the density dependence of the asymmetry term. The data are also shown as the solid points. 
The various hues in the figure are described in the text. The excitation energies for the 
fragments are taken directly from the SMF calculations. Lower panel: The data are the 
same as in the middle panels. The curves are the calculations obtained when the excitation 
energies of the primary fragments are reduced by a factor of two. 

Figure 7: Dependence of the scaled function S(N) on the density dependence of the asym- 
metry term. The left panel provides a comparison between values for S(N) computed from 
the data (solid line) and the calculated primary (points about dashed line) and final (points 
about dashed-dotted line) distributions obtained for the super stiff asymmetry term. The 
right panel provides a comparison between values for S(N) computed from the data and the 
calculated primary and final distributions obtained for the soft asymmetry term. 
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